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Recent high precision measurements of the CMB anisotropies performed by the BOOMERanG and 
MAXIMA- 1 experiments provide an unmatched set of data allowing to probe different cosmological 
models. Among these scenarios, motivated by the recent measurements of the luminosity distance 
versus redshift relation for type la supernovae, is the quintessence hypothesis. It consists in assuming 
that the acceleration of the Universe is due to a scalar field whose final evolution is insensitive to 
the initial conditions. Within this framework we investigate the cosmological perturbations for 
two well-motivated potentials: the Ratra-Peebles and the SUGRA tracking potentials. We show 
that the solutions of the perturbed equations possess an attractor and that, as a consequence, 
the insensitivity to the initial conditions is preserved at the perturbed level. Then, we study the 
predictions of these two models for structure formation and CMB anisotropies and investigate the 
general features of the multipole moments in the presence of quintessence. We also compare the 
CMB multipoles calculated with the help of a full Boltzmann code with the BOOMERanG and 
MAXIMA- 1 data. We pay special attention to the location of the second peak and demonstrate 
that it significantly differs from the location obtained in the cosmological constant case. Finally, we 
argue that the SUGRA potential is compatible with all the recent data with a standard values of 
the cosmological parameters. In particular, it fits the MAXIMA-1 data better than a cosmological 
constant or the Ratra-Peebles potential. 



I. INTRODUCTION 

Recent measurements of the luminosity distance versus 
redshift relation for type la supernovae [jjj, if confirmed, 
are compatible with an expanding (accelerating) universe 
driven by a new type of matter whose equation of state 
p = ujp is characterized by a negative u>. One of the pos- 
sible explanations is the existence of a non-zero vacuum 
energy, i.e. a "cosmological constant". Another prag- 
matic possibility which has been proposed is to assume 
the existence of a yet unknown mechanism guarantee- 
ing that the true cosmological constant vanishes, the re- 
maining energy density being then due to the presence 
of a scalar field, the quintessence field, almost decoupled 
from ordinary matter [Q-|5| . The main difference between 
a quintessence fluid and a cosmological constant comes 
from their equation of state where wa = — 1 for a cosmo- 
logical constant and — 1 < ujq < for the quintessence 
fluid. 

One of the puzzles in the interpretation of these data 
is the extremely small value of the energy density due 
to the new form of matter. From the point of view of 



particle physics a vanishing value for the cosmological 
constant is one of the major challenges ||. At present 
there is no known mechanism which prevents the vacuum 
energy from picking large values due to radiative cor- 
rections and one expects typically a contribution equal 
to {he 1 2) f dk/c/(27r) 3 ~ hck^ ax /{16ir 2 ), where /c max is 
a cut-off which can naturally be taken as the Planck 
wavenumber. This gives a contribution which is 120 or- 
ders of magnitude above the observed one. One possibil- 
ity which is often advocated is the presence of some global 
supersymmetry (SUSY) which would guarantee that the 
energy of the vacuum is zero. Unfortunately SUSY has to 
be broken to take into account the absence of experimen- 
tal evidence in favour of particle superpartners leading 
to a natural contribution to the vacuum energy of order 
Afg USY where MgusY is the SUSY breaking scale esti- 
mated around 1 TeV [Q . The measurement of a vacuum 
energy some 60 orders of magnitude below this expected 
value indicates that some new physics must be at play 
here. 

In the quintessence hypothesis, the small vacuum en- 
ergy density is due to the rolling down of the quintessence 
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field Q along a decreasing potential. A typical poten- 
tial is the Ratra-Peebles potential V(Q) = A 4+a /Q a [§. 
From the particle physics point of view one would like 
to justify the existence of the quintessence field. Sev- 
eral natural candidates have been ruled out such as the 
axion-dilaton field || , the moduli fields of toroidal com- 
pactifications in string theory || and finally the meson 
fields of supersymmetric gauge theories |L(J. Neverthe- 
less, it seems reasonnable to expect that SUSY will play 
a role in the solution. Within this framework it is a 
matter of fact that the quintessence field must be part 
of supergravity (SUGRA) models This comes 

from the large value Q ~ mpi of the field at small red- 
shift which implies that SUGRA corrections cannot be 
neglected. In 11 1 an effective theory approach has been 
used to deduce the general form of quintessence SUGRA 
potentials, they are of the type 
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where k = 8ttG, G being the Newton constant, and where 
the exponential factor comprises the SUGRA corrections. 
A and a are free parameters. The fine-tuning is not too 
severe as for typical values a = 6 the scale A ~ 10 6 GeV 
is compatible with high energy scales. Notice that the 
SUGRA corrections become relevant towards the end of 
the evolution and decouple at small Q <C mp\. Different 
types of potentials can be distinguished because they lead 
to different values of the equation of state parameter. For 
example, for a = 11, the Ratra-Peebles potential is such 
that wq ~ —0.29 whereas the SUGRA potential gives 
cu Q ~ -0.82 @ (for fiq = 0.7). 

It is also worth noticing that there exists quintessence 
models where the field is non-minimally coupled with the 
metric. Such models induce time-variation of the New- 
ton constant and are therefore already constrained, for 
example by observations in the solar system or by pul- 
sar timing measurements p2| , p^| . They lead to the same 
tracking behaviour, as stressed in Refs. [jl4|,[l5|, as soon as 
the coupling term is proportional to a power of the poten- 
tial. However, some important differences occur when the 
field starts dominating; for example its effective equation 
of state can reach extreme values such that u> ~ —3 Jl6| . 
Also, these models can lead (especially in the context of 
quintessential inflation |l7j]) to clear observable features 
in the gravitational waves spectrum jl8| . 

In view of the numerous phenomenological successes 
of quintessence, it is relevant to deduce its consequences 
for Cosmic Microwave Background (CMB) anisotropics 
and structure formation. The aim is two-fold. First, we 
have to study whether quintessence leads to acceptable 
scenarios and, second, we have to learn how we could use 
high-precision measurements recently obtained by the 
BOOMERanG @-|l| and MAXIMA-1 HQ experi- 
ments or to be performed in the near future by NASA's 
Microwave Anisotropy Probe (MAP) satellite Q, ESA's 
Planck satellite pa] or the Sloan Digital Sky Survey 



(SDSS) |2(| to put constraints on the quantities char- 
acterizing quintessence like or ojq. The second possi- 
bility has of course already been investigated for the cos- 
mological constant case. For example, the fraction Q\ of 
the critical density is not determined entirely from the 
supernovae data. Indeed, the data from the supernovae 
observations are degenerate in the plane (f2 m , £Ia), where 
f2 m is the matter (i.e. cold dark matter plus baryons) 
component preventing a clear cut determination of the 
fraction The situation changes drastically if one in- 
cludes the measurements of the CMB anisotropics |27| ] 
(even without the BOOMERanG or MAXIMA-1 data). 
In that case, the degeneracy is removed leading to a prob- 
able 70% of the total energy density of the universe car- 
ried by the negative pressure fluid while the remaining 
30% are the matter components ensuring that Qq = 1 in 
agreement with a spatially flat universe. This conclusion 
can be drawn from the measurements of the location of 
the first Doppler peak. This result has been confirmed 
by other measurements p8|-|30|. Another use of com- 
bined data has been to put constraints on the equation 
of state parameter. However, this has been done only for 
constant or for very simple time-dependent ujq ]3"H-|33]| . 

CMB anisotropics and the power spectrum are calcu- 
lated with the help of the theory of cosmological per- 
turbations. Cosmological perturbations in the presence 
of quintessence have been studied by Ratra and Peebles 
but only in the tracking regime || . CMB multipoles mo- 
ments and/or the power spectrum have already been cal- 
culated for the Ratra-Peebles potential in Ref. J 34 and 
for other models of quintessence in Refs. [p5|-^9|. One 
important issue is to understand whether the final evo- 
lution of the various perturbed quantities depend on the 
initial conditions imposed at reheating (of the inflation- 
ary type or not). Another way to put the same problem 
is the following: do the multipole moments depend on the 
value of 5Q(r]i) and 5Q'(?]i) at initial time? In Ref. [j34| , 
it was noticed that the answer to this question is no but 
no explanations were provided. Here, we confirm the re- 
mark of Ref. |34| and show that this is due to the fact 
that the perturbed Einstein equations also possess an at- 
tractor which renders the multipole moments insensitive 
to the initial conditions. 

One of the main purpose of this article is the study of 
the general properties of the multipoles moments of the 
CMB anisotropies in presence of the quintessence field. 
We present the CMB multipole moments for the Ratra- 
Peebles potential and, for the first time, for the SUGRA 
tracking potential. In addition, we also display the mat- 
ter power spectrum for these two models. Recently, it 
has been shown by Kamionkowski and Buchalter ^] 
that the location of the second peak in the CMB power 
spectrum is an efficient way of revealing some features 
of the dark energy sector. Therefore, we pay special at- 
tention to this question. In particular, in Ref. p0[ , only 
the cosmological constant case was studied and it was 
argued that the quintessence case (the authors refer to 
the Ratra-Peebles potential) must not differ significantly 
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from the cosmological constant case. In the present ar- 
ticle, we demonstrate that this is not the case and that, 
as a matter of fact, quintessence leads to a different loca- 
tion (denoted, in the following, I2) of the second peak. In 
addition, we show that the location of the second peak 
in the quintessence case and in the cosmological constant 
case can be easily distinguished. Following Ref. jlO|], wc 
display the contour plots of I2 in the plane (£l m ,h) for 
the Ratra-Peebles and SUGRA tracking potentials. 

The article is organized as follows. In section |l[ we 
give a description of the background evolution in terms of 
two physical quantities: the equation of state parameter 
ojq and the sound velocity c 2 q. In section III, we study 
the cosmological perturbations for the quintessence field. 
In section [V , we present the results of full numerical cal- 
culations with the help of a Boltzmann code developped 
by one of us (A.R.) for the CMB anisotropics and power 
spectra in the case of the Ratra-Peebles and SUGRA 
potentials. Then, detailed comparisons with the recent 
BOOMERanG @|IJ and MAXIMA- 1 g|| data are 
performed. We end with our main conclusions in sec- 
tion |v|. 



II. THE BACKGROUND EVOLUTION 

We suppose that the Universe can be described by a 
Friedman-Lemaitre-Robertson- Walker metric the space- 
like sections of which are flat 



oj' q = -3H{1 + w Q )(c s 2 Q - u Q ), 



(4) 



ds 2 = a 2 (?7)(-d?7 2 + Sijdx'dx*). 
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In this equation, 77 is the conformal time related to the 
cosmic time by a(rj)dr] = dt. The matter content is as fol- 
lows. The Universe is filled with a mixture of five fluids: 
photons (7), neutrinos (is), baryons (b), cold dark mat- 
ter (cdm) and a scalar field Q named quintessence. The 
stress energy tensor of each of these species is the one of 
a perfect fluid, = (p + p)u^u v +pg^, where it M is the 
4-velocity of the fluid. The energy density and the pres- 
sure of the scalar field are given by pq — \{Q' /a) 2 +V(Q) 
and pq = i(Q'/a) 2 — V(Q), where V(Q) is the poten- 
tial of quintessence whose shape will be very important 
in what follows. Each fluid is also characterized by its 
equation of state pi = LO\p\ where i = 7, v, b, cdm or 
Q. We have oj 1 = oj v = 1/3 and oj^ — u c dm = 0. 
The case of ojq is more complicated since this is a time- 
dependent function such that — 1 < ujq < +1. Its ex- 
pression reads ojq = 1 — 2V(Q)/pQ. The fact that ojq is 
a time-dependent function directly comes from the fact 
that, for a scalar field, the sound velocity defined as Q] 



4a 2 dV{Q) 
3HQ' dQ 



K 2 w +i )' (3) 



is not equal to the equation of state parameter ojq. As a 
consequence ojq has to change in time as revealed by the 
following equation 



unless ujq = — 1. 

The evolution of the Universe can be calculated with 
the help of the Friedman and conservation equations 



-o At — —j- > pi, 



(5) 



Pi 



-3H(l + 0Ji)pi, i = 7, v, b,cdm or Q, (6) 



where mp\ is the Planck mass and TL = a' /a is related 
to the Hubble constant by the equation H = H/a. The 
equations of conservation simply express the fact that the 
energy is conserved for each species which do not interact. 
The equation of conservation of the quintessence field can 
also be written as the Klein-Gordon equation 



Q" + 2HQ' + a 2 ^- = 0. 

dQ 



(7) 



We now need to give the last piece of information neces- 
sary to have a complete description of the system, i.e. the 
shape of the potential V(Q). In order to be an interest- 
ing theory and to represent an improvement over the cur- 
rent situation, quintessence has to address the following 
four problems: the fine-tuning problem, the coincidence 
problem, the equation of state problem and the model 
building problem. The fine-tuning problem amounts to 
understanding whether one can have Qq ~ 0.7 with the 
free parameters of the potential taking "natural" values, 
i.e. close to the energy scale of the theory under consider- 
ation. The coincidence problem is the question of the ini- 
tial conditions: does the final value of pq strongly depend 
on the chosen initial values of Q and Q'l The equation 
of state problem is the question of the value of ojq. In or- 
der to be compatible with observational data, it should 
be such that —1 < ojq < 0. According to recent pa- 
pers, even more stringent restrictions can be put, namely 
-1 < ojq < -0.6 H or evcn -1 < < -0.8 || . In 
particular, this already rules out a network of cosmic 
strings since the corresponding fluid has an equation of 
state parameter equal to —1/3. Finally, the model build- 
ing problem consists in justifying the shape of the poten- 
tial from the high energy physics point of view. Many 
different shapes of potential which allow, at least par- 
tially, to solve these problems have been investigated in 
the litterature and Table | summarizes these proposals. 
In particular, the first possibility has been studied thor- 
oughly in the past years. In this article, we will mainly 
concentrate on the Ratra-Peebles potential [Q and the 
SUGRA tracking potential 0Jl|]. 

Let us briefly see how the four questions evoked previ- 
ously can be addressed with these potentials. 



A. The fine-tuning problem 

Let us start with the fine-tuning problem which is 
clearly a delicate question. This problem is crucial M 
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TABLE I. Quintessence potentials that have been used in 
the litterature. 



for the cosmological constant. Indeed, from very sim- 
ple high energy physics considerations, one typically ex- 
pects pA — rripY ~ 10 76 GeV 4 whereas one measures p\ ~ 
fiAPc — 10~ 47 GeV 4 since the critical energy density is 
p c ~ 8.lh 2 x 10- 47 GeV 4 . Do we gain something in the 
case of quintessence? This question is controversial. For 
example in Ref. p2[ , the authors clearly answer no and 
write "Two proposals to explain these observations are 
a non-vanishing cosmological constant or a very slowly 
rolling scalar field, often dubbed quintessence. Both pro- 
posals, however, are plagued with formidable fine tun- 
ing problems." However, one should look more carefully 
at this point. To illustrate this issue, let us consider 
the general argument given against quintessence. If we 
consider the potential V(Q) — (m 2 /2)Q 2 then the mass 
of such a field, which is also the only free parameter of 
the potential, should be m = ^2flQp c /mpi ~ 10 -33 eV, 
a very tiny mass indeed. Justifying such a value for 
the free parameter m is probably the same problem as 
justifying a very low value for p\. However, such a 
model has never been advocated for the quintessence 
field. As already mentioned above, one typically con- 
siders models such that V(Q) — A 4+a /Q a . This changes 
the argument. Now, the free parameter of the theory 
is A. In order to have pq = CIqp c today, one has 
A ~ 10 11 GeV, for a = 11. This time, the free parame- 
ter of the theory has a value comparable to the natural 
scales of high energy physics. Therefore, something has 
been gained and it seems unfair not to emphasize this 
point. On the other hand the mass of the field is given 
bym= a(a + l)ClQp c /m Pi ~ 10 _33 eV but this number 
should be interpreted completely differently. Here the 
mass m is just a "by-product" and its value is naturally 
very small without any artificial fine-tuning of A. Of 
course the very small value of the mass implies that the 
quintessence field is almost completely decoupled from 



the other matter fields. This renders the model building 
issue even more acute. 



B. The coincidence problem 

The coincidence problem as formulated in the intro- 
duction, i.e. the dependence upon the initial conditions, 
is solved because the Klein-Gordon equation possesses 
an attractor. In order to prove this property, we have 
to rely either on numerical calculations or on approxi- 
mate methods. All the plots and numerical estimates 
displayed in this article will be made with the help of 
numerical calculations. However, it is always useful to 
understand the tracking property by means of analytical 
methods and we now turn to this question. It is conve- 
nient, for analytical calculations, to consider that there 
is in fact only one "background" fluid with a time de- 
pendent equation of state such that wb = 1/3 during the 
radiation dominated epoch and wb = during the matter 
dominated era. In addition to the background fluid, we 
assume that there also exists the quintessence scalar field 
field. Following the treatment of Ratra and Peebles 0, 
it will be considered that this scalar field is a test field. 
This is a good approximation since this field must be sub- 
dominant in particular during Big Bang Nucleosynthesis 
(BBN) in order not to modify the behaviour of the scale 
factor and, as a consequence, not to spoil the success of 
BBN. This means that the behaviour of the scale factor is 
essentially determined by the background fluid and that 
pi ~ pb- This hypothesis breaks down at very small 
redshift when quintessence starts dominating the matter 
content of the Universe. Since quintessence is only a test 
field which does not interact with the background fluid, 
the scale factor and the quantity Ti. can be written as 

a( V ) oc r, 2 /( 1+3 ^) , n{n)= 2 . (8) 

For the sake of illustration, let us now consider the radia- 
tion dominated era where o>b = 1/3. Under the previous 
assumptions, the Klein-Gordon equation has a particular 
solution given by 

Q P = g r/ 4/(Q+2) , (9) 

where Qq is a constant which depends on the free pa- 
rameters of the potential, i.e. A and a. The tracking 
behaviour is revealed by the behaviour of small pertur- 
bations around Q p . Let us introduce the new time r 
defined by rj = e T and define u and p by Q = Q p u and 
p = du/dr. The Klein-Gordon equation, viewed as a 
dynamical system in the plane (p, u), possesses a critical 
point (0, 1) and small perturbations around this point 
5u, 5p obey the following equation 
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FIG. 1. Evolution of the kinetic energy (red line) and of the 
potential energy V(Q) (green line) from z\ = 10 30 to z = 0. 



Solutions to the equation det(M — XI) = 0, where M is 
the matrix defined above, are given by 



2(a + 2) 2(a + 2) 



\/l5a 2 + 108a + 92. (11) 



The real part of X± is always negative and the critical 
point is a spiral point. Therefore, every solution will 
tend to Q = Q p after an intermediate regime: Q = Q p is 
an attractor and no fine-tuning of the initial conditions 
is required. 

Before reaching the attractor, the quintessence field 
undergoes different regimes that we are now going to de- 
scribe. These regimes are in fact characterized by two 
physical quantities already introduced previously: the 
equation of state parameter uiq and the sound velocity 
CgQ. We study the case of an "overshoot", in the termi- 
nology of Ref . H , since this corresponds to initial condi- 
tions that are physically more relevant (in particular this 
includes the case of equipartition, i.e. pq ~ 10 ini- 
tially). We also assume that the background is radiation 
dominated ub = 1/3. 

Initially, the kinetic energy dominates the potential en- 
ergy, i.e. Q' 2 /(2a 2 ) 3> V(Q). This means that the energy 
density redshifts as pq cx 1 /a 6 and that the equation of 
state parameter is ujq = 1. As a consequence, due to 
the constancy of luq and Eq. (Q) (and also ujq ^ —1), we 
have c 2 q = 1 as well. The scalar field itself evolves like 



Q = Qi 



A 
a 



(12) 



where Qf and A are constant. These constants are such 
that the term A/a becomes rapidly small in comparison 
with the frozen value Qf and we have the amusing situ- 
ation that the field can be (almost) considered as frozen 
even if the kinetic energy still dominates. This is illus- 
trated in Fig. [l]. 

As a consequence, during this regime the potential en- 
ergy is also almost constant except at the very beginning. 



Using the definition of ujq and c 2 q, see Eq. (||), we deduce 
that, during the kinetic regime, we have 



ujq — 1 cx a 



1 cx a° 



(13) 



The fact that, in the parametrization adopted here, the 
scale factor is very small during the kinetic regime ex- 
plains that there is no contradiction between these equa- 
tions and the values of ujq and c 2 q deduced above. 

Since the kinetic energy decreases while the potential 
energy is almost constant, the kinetic regime cannot last 
forever. When the potential energy becomes larger than 
the kinetic one, the equation of state parameter suddenly 
jumps from +1 to —1 while the sound velocity still re- 
mains equal to +1 since Eq. (|]) does not imply a change 
of this quantity in the case loq — — 1. The fact that the 
equation of state parameter changes before the sound ve- 
locity is explained by Eq. (p"3|). We call this regime the 
transition regime. During this regime, the kinetic en- 
ergy still redshifts as 1/a 6 and V(Q) is approximately 
constant but of course now Pq — V(Q). 

Due to the second of Eq. ( |l3| ) , the sound velocity has 
also to change at some later time. This implies that 
the quintessence field can no longer behave according 
to Eq. (|l2[). This is the starting point of the potential 
regime. In order to study the behaviour of the system 
in this regime, we need to find an expression for the sec- 
ond derivative of the potential. Differentiating once the 
definition of the sound velocity, Eq. (0), we arrive at 



d 2 V(Q) 
dQ 2 



?W^c 2 ' 



+ (c s 2 



H 



1 



n 2 



(3c s 2 Q + 5) 



(14) 



No approximation has been made in the derivation of this 
relation. This formula generalizes Eq. (3) of Ref. This 
formula will turn out to be very useful when we study the 
perturbations in the next section. With the scale factor- 
given by Eqns. (jq), this relation can be re- written as 



2 d 2 V(Q) 
3H 2 dQ 2 



lc 2 
H 



3(c 2 Q ^l)( WB + C , 2 Q + 2) 



(15) 



In the regime we are interested in, the r.h.s of the previ- 
ous formula is small. The only way to satisfy this relation 
is to ensure that the sound velocity changes to the con- 
stant c 2 q = — 2 — ujb- This gives c 2 q = —7/3 for the 
radiation dominated era. This evolution is displayed in 
Fig. § 

The fact that the sound velocity is a constant implies 
that the factor (4a 2 )/ (3HQ')dV(Q)/dQ is also a con- 
stant. Therefore, the behaviour of the quintessence field 
is now given by 

Q = Q f + Ba 4 , (16) 
which implies that the kinetic energy redshifts as a 4 . 
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FIG. 2. Evolution of the equation of state parameter loq 
(green line) and of the sound velocity c 2 q (red line) from 
z< = 10 30 to z = 0. 
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1 
1 
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TABLE II. Summary of the different regimes described in 
the text. 



Again this regime cannot last forever since the kinetic 
energy increases while the potential energy still remains 
constant. At some later time, both contributions become 
equal and ujq and c^q have to change once more. This is 
the end of the potential regime and the beginning of the 
tracking regime which has already been described above. 
The quantities pq, pq, V and the kinetic energy reach a 
fixed ratio such that 



2 — auJB 



(17) 



The definitions of the different regimes and the corre- 
sponding evolutions of the physically relevant quantities 
are summed up in Table [l| 



C. The equation of state parameter problem 

The third question evoked previously was the question 
of the value of the parameter ujq today. As already men- 
tioned, this is an important issue since constraints on 
this quantity are already available. This problem is also 
solved by quintessence in the sense that we always have 
— 1 < ujq < 0. Here, however, it is relevant to distinguish 
between the Ratra-Peebles potential and the SUGRA po- 
tential. With the first potential, it seems difficult to reach 
sufficiently small value of ujq. On the other hand, this 
is automatically achieved in the second case. The reason 



for this is the presence of the factor exp(ft<3 2 /2) in the 
potential, a generic feature of SUGRA-based potentials, 
which drives ujq towards —1. For a = 11 and fl c d m — 0.3, 
the prediction is ujq ~ —0.82 a value in agreement with 
the current data ]Tc| , pd| . 



D. The model building problem 

From the particle physics point of view, one would 
like to justify the existence of the quintessence field and 
the shapes of the (so far) phenomenological potentials. 
Several attempts have already been made in the frame- 
work of supersymmetric field theory. In particular, it was 
shown by Binetruy ||] that the Ratra-Peebles potential 
can be recovered in the context of global SUSY. How- 
ever, as already mentioned, SUGRA corrections must be 
taken into account and this implies that the correspond- 
ing potential can be of the type of the SUGRA tracking 
potential displayed in Eq. (|l]) which leads to a better 
agreement with the available data. 

Nevertheless, it should be clear that considerable prob- 
lems remain to be addressed in order to reach a satisfac- 
tory situation . Maybe the most crucial question is 
the problem of SUSY breaking. SUSY must certainly be 
broken but the models evoked previously do not take into 
account this basic fact. This could have dramatic conse- 
quences and modify the shape of the potential which is so 
important in order to solve the three previous problems. 



III. THE COSMOLOGICAL PERTURBATIONS 

We now turn to the study of the cosmological pertur- 
bations. A detailed study has already been performed 
by Ratra and Peebles in Ref. (|] but only for the track- 
ing regime. Cosmological perturbations in a fluid with a 
constant negative equation of state parameter have been 
investigated in Ref. (IjJ. In this article, we study the 
cosmological perturbations (in the long wavelength ap- 
proximation) in all the regimes previously described and 
point out some additional properties. The evolution of 
the cosmological perturbations mainly depends on the 
equation of state parameter and the sound velocity. We 
have shown in the previous section that they can be con- 
sidered as constant in each regime. This will simplify the 
analysis a lot. 

The fate of the perturbations depends on the initial 
conditions. It has been noticed for the first time in 
Ref. H] that "the observable fluctuation spectrum is in- 
sensitive to a broad range of initial conditions, including 
the case in which the amplitudes of SQ, SQ' are set by in- 
flation" . In that paper, the authors choose SQ = SQ' = 
initially (in the synchronous gauge). We demonstrate, in 
this section, that the insensitivity of the spectrum de- 
scribed in Ref. p4| has an origin similar to the insensi- 
tivity of the background properties with respect to the 
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initial conditions Q and Q' , namely the presence of an at- 
tractor for the perturbed quantities. We prove that dur- 
ing all the four regimes undergone by the quintessence 
field, the attractor is characterized by a "spiral fixed 
point" as it is the case for the background. 



A. General framework 

Without loss of generality, the perturbed line element 
can be written in the synchronous gauge. In this class of 
coordinates systems, scalar perturbations are completely 
described by two arbitrary functions. The spatial depen- 
dence of the perturbations is given by X(x > ) which is the 
eigenfunction of the Laplace operator on the flat space- 
like hypersurfaces. There exists two ways to construct a 
two rank tensor from a scalar function : either by mul- 
tiplying it by the spatial background flat metric 6ij or 
by differentiating it twice. The two arbitrary functions 
mentioned above are simply the coefficients of these two 
terms in a Fourier expansion. Therefore, the perturbed 
metric can be expressed as [pOl 



ds 2 = a z {ri)\ -dry 



I + h{rj)X \ Sij 

dzW'l. (18) 



In this equation, the dimensionless quantity k is the co- 
moving wavevector related to the physical wavevector 
k phys through the relation k phys = \s./a{rj). As a conse- 
quence of Einstein equations, perturbations in the metric 
are coupled to perturbations in the different matter com- 
ponents. We choose to write the perturbed stress-energy 
tensor according to fKJ 
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er 

1*3 = 



T l o = -KX 



(19) 
(20) 



where we have assumed that the longitudinal pressure 
pi vanishes for each component. As for the background, 
one considers that the Universe is filled with two fluids: 
the background fluid, an hydrodynamical perfect fluid 
which is either radiation or dust (again, the correspond- 
ing quantities will carry the index B) and a scalar field 
Q describing the quintessence field (in this case the cor- 
responding quantities will carry the index Q). The per- 
turbed Einstein equations which govern the evolution of 
the quantities h and hi are given by: 



mti 



ti, 



\- k 2 h — Tih[ — k£ib + kciq, 

ti = <B + <Q, 

- h" - 2Hh' = npiB + npiQ, 

mti, - k 2 h = o. 



(21) 
(22) 
(23) 
(24) 



Finally, it turns out to be more convenient to work with 
the density contrast 8 and the velocity divergence 6 de- 
fined by the equations: 



ei 
a 2 e 



(25) 



In the following, we study analytically the time evolution 
of the density contrast for the background fluid and for 
quintessence in the long wavelength limit. 



B. The background fluid 

The equations satisfied by the background density con- 
trast and divergence can be obtained either from com- 
binations of the Einstein equations (pl[-p4|) or, more 
directly, from the conservation of the perturbed back- 
ground fluid stress-energy tensor (since the background 
fluid and quintessence only interact gravitationally) . 
They read: 

<5 B + a(l + u^)0 B + ^^(3/i' -h'i) = 0, (26) 



6' B + (2 - 3w b )H0b - 7Y 



(1 + uj B )a 



sB S B = 0. 



(27) 



These two equations are equivalent to Eqns. (7.15) 
and (7.16) of Ref. |). From them, we can derive the 
relation 



-2 



9k 2 r 2 



3h" - ti{ = 8'i + 2(1 - 3uj B )a'8 B - S -^S B , 

1 + <jj b 1 + oj B 

(28) 

where we have assumed that oj b is a constant. On the 
other hand, from the Einstein equations we get 



-(3h" - h'l) - H(3ti - ti t ) 



3W 



(1 + 3c} B )Sl B 5 B + (1 + 3c, 2 Q )fi Q <S Q 



(29) 



where c 2 q = piq/eiq which needs not to coincide with 
the definition of c 2 q. In order to derive the formula satis- 
fied by the density contrast of the background fluid in the 
long wavelength limit, we neglect the term proportional 
to k 2 in Eq. ( p8| ) and we use the fact that Qq -C S!b- 
Then, straightforward manipulations lead to 



6& + H6' B - -W 2 (1 + 3^ b )(1+^b)£b 
= -3H^b(1 + u B )aQ B , 



(30) 



where we used the fact that cj B = uj b for an hydrody- 
namical fluid. This equation shows that the evolution of 
the background density contrast is essentially unaffected 
by the presence of quintessence. This is of course an 
expected result since we have assumed Qq <C tt B . The 
general solution to Eq. (|30|) can be easily found and reads 
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TABLE III. Time dependence of the background fluid den- 
sity contrast during radiation and matter dominated era 



a \ A I a 

+ A 2 [ — 



Sb(v) = M \ — 
a 



«o 



w B (l + w B )(l + 3LU B )a a 9 B0 r) ( a \ (9 " B 1)/2 



(l-w B )(l + 6w B ) 
where we have defined 
(1 - 3lu b ) 



rt 



(31) 



1 



±- V(l - 3w B ) 2 + 24(1 + w B )(l - 3wb). 



(32) 



The results for the radiation dominated and matter dom- 



inated epochs are summarized in Table III . These results 
are consistent with those obtained in Ref. |5^]. In par- 
ticular, it can be shown that the branch <5 B oc a x ~ corre- 
sponds in fact to a residual gauge mode, i.e. there exists 
a synchronous system of coordinates such that this mode 
can be removed and therefore must not be considered as 
a physical mode. 



C. Quintessential perturbations 

We now describe how the long wavelength quintessen- 
tial perturbations evolve with time. A similar study has 
already been performed by Ratra and Peebles but only on 
the tracking solution. We give here a complete descrip- 
tion of the evolution of the quintessence density contrast 
in the four regimes defined in the previous section. In 
addition, we prove that there exists an attractor for the 
perturbations as it is the case for the background solu- 
tion. As a consequence, the final value of the density 
contrast is always the same whatever the initial condi- 
tions are. 

In order to obtain the fundamental equations to be 
solved, we can proceed as for the background fluid. How- 
ever, it is important to notice that the link between the 
perturbed energy density and the perturbed pression, 
which is just a constant for the background fluid, is more 
complicated in the case of quintessence. In general, we 
can write piq = c^eiQ + ci 2 t5S where the second term 
represents entropy perturbations. In the synchronous 
gauge, we obtain 



Piq 



C sQ £ 1Q 



1 



+ (1-<* Q )-W + Hh\) 



(33) 



We can now establish the equations satisfied by the 
quintessence density contrast and divergence. The con- 
servation of the perturbed stress-energy tensor leads to 



S' Q + 3a(c 2 Q - u Q )S Q + a(l + u Q )0 Q 



+ ±(1 + w Q )(3fc' - h[) = ^-(c s 2 Q - l){h'{ + Hh[), (34) 
9' Q + (2 - 3^ Q ) 



s Q 



[l + WQ)a 
(l-c 2 Q )fc 2 



1 + cjQ w 3aH 2 (l+uJ Q )n Q 



{h'{+Hh[). (35) 



In these two equations, no approximation have been 
made: they are valid for any wavenumber, any equa- 
tion of state parameter and sound velocity. In practice, 
it turns out to be more convenient to use the perturbed 
Klein-Gordon equation to analyse the problem. This can 
be obtained directly from the first of the two previous 
equations if one notices that the quantities describing 
the perturbed scalar field stress-energy tensor can be ex- 
pressed in terms of the perturbed scalar field SQ(r],x) 
according to 



e iQ 

£q 

PlQ = Q' 



-Q'SQ, 



} dV(Q) 
' dQ • 

? dV(Q) 
' dQ • 



(36) 
(37) 
(38) 



Inserting the corresponding expression for the density 
contrast and the divergence in Eq. (p54), we get 



SQ" + 2H5Q' + 



1 dQ^ 



SQ+ y(3/i'-/i0 = 0. 

(39) 



This is similar to Eq. (7.20) of Ref. §. One can check 
that Eq. ( |35| ) is automatically verified since it is equiv- 
alent to the unperturbed Klein-Gordon equation (times 
an unimportant factor). Using Eq. (26) to express the 
factor 3h' — h', and neglecting the k 2 term, we arrive at 



SQ" + 2H5Q' 



2 d 2 V(Q) 
1 dQ 2 



5Q = Q'a8 B 



Q' 



1 + U B 



(40) 



We are now going to analyse this equation in detail. We 
now need to utilize the general expression for the second 
derivative of the potential, Eq. (]14|). On the tracking 
solution, we have loq = c 2 q andwQ = (— 2+ao; B )/(2+a) 
and this equation reduces to 



d 2 V{Q) 
dQ 2 



9 g2 a + 
2 a 



1 



-(1 



(41) 



For our purpose, as proven in the previous section, it 
is sufficient to consider a regime where c 2 q is constant 



and where the scalar field is a test field, 
conditions, we obtain 



Under these 
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d 2 V(Q) 
dQ 2 



-ff 2 (l- C 2 )(6 + 3c, B + 3 C 2 ). (42) 



Let us now concentrate on the homogeneous part of 
Eq. (p0|). Using the previous equation, it can be ex- 
pressed as 



SQ 1 ' 



3ui 
1 



(1 



3w B ) 2 V 2 



(l-c 2 Q )(6 



3w B + 3c s 2 Q )5Q = 0. (43) 



This linear equation can easily be solved: its solutions 
are just power law of the conformal time. However, in 
order to show explicitly the complete analogy with the 
background attractor, we choose to analyse it in a rather 
roundabout way. Let us proceed exactly as for the unper- 
turbed Klein-Gordon equation [see the discussion around 
Eq. (|lC)|)]. We define the time r by r\ = e T and intro- 
duce the quantity Su and Sp defined by Su = SQ and 
Sp = d(SQ)/dT. Then, Eq. (pj|) can re-expressed as 



d_fSp 
dr V Su 



I 3(c^ B -l) 
= l+3w B 

\ 1 



S£f(2 + -b + c s 2 q ) 





(44) 



The form of this equation clearly shows the complete 
analogy with Eq. ( |l0| ) . The eigenvalues of the system are 
found by solving the equation det(M — XI) = 0, where 
M is the matrix defined above and / the identity matrix. 
Straigthforwards calculations show that the solutions are 
given by 



A+ = - 



3 wb - 1 



2 1 + 3lj b 



1± Wl+4 



1 



(WB - 1) ! 



r (2 + ^ B +c 2 , 



(45) 



Of course, this is just a simple rephrasing of the fact that 
the solution of Eq. (||) is SQ oc A + t] X + + A-i] X - . The 
presence of an attractor is linked to the negative sign 
of the real part of \±. It is easy to see that the real 
part is always negative in all four regimes, in particu- 
lar this is true for any value of a. This is displayed in 
Fig. U in the plane (wb,c 2 q). The green and purple re- 
gions are the regions where these real parts are negative. 
The green region is the region where the argument of the 
square root is negative, i.e. where the square root is an 
imaginary number. The exact "trajectories" of the sys- 
tem for the usual tracking potential (short line) and for 
the SUGRA tracking potential (long line) are also shown 
for the case a = 11. They have been obtained by full 
numerical integration. The remarkable property is that 
these trajectories are always in the stable region. This 
means that, in each region, the system tends to an at- 
tractor which is given by the inhomogencous part of the 
perturbed Klein-Gordon equation. The system starts at 
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FIG. 3. Stability region for the quintessential perturba- 
tions. The green (lightgrey) and the purple regions (dark- 
grey) constitute the stable region. The green region is the 
region where the square root is an imaginary number. Two 
trajectories of the system for the tracking potential (shortest 
line) and the SUGRA tracking potential (longest line) are also 
displayed for the value a = 11. 



uj b = 1/3 and goes from c 2 q = 1 to c 2 q = —7/3. Then, 
the system approaches the transition to the matter dom- 
inated era and leaves the vertical line. Finally, it stops 
when the redshift vanishes at wb — —0.29 for the track- 
ing potential and at uo-q ~ —0.82 for the SUGRA tracking 
potential. The two lines separate when the exponential 
factor becomes important in the SUGRA tracking poten- 
tial. 

The conclusion is that the final value of the 
quintessence perturbations is insensitive to the initial 
conditions, a property completely similar to what has 
been shown in Ref. J5j for the background. Strictly speak- 
ing, this property has been demonstrated for long wave- 
length modes only. However, we have checked by numer- 
ical calculations that this is also true for shorter wave- 
length modes. Having proven that the final result does 
not depend on the initial conditions of the quintessence 
perturbations, we can now proceed further and embark 
in a rather detailed study of the CMB anisotropies pre- 
dictions in the presence of quintessence. 



IV. PREDICTIONS FOR THE POWER 
SPECTRUM AND THE MULTIPOLE MOMENTS 

The presence of cosmological perturbations induces di- 
rectional variations in the CMB photon redshift. This is 
the so-called Sachs- Wolfe effect |52|. Since these varia- 
tions are the same regardless of the wavelength of the 
photons, they translate into variations in the tempera- 
ture of the black body on the celestial sphere. Their 
amplitude has been measured by the COBE satellite and 
is of the order of magnitude <5T/T ~ 10" 5 [|§. The de- 
tailed angular structure of the CMB anisotropies is usu- 







ally characterized by the two-point correlation function 
which can be expanded according to 

^( ei )^(e 2 )\ = i-£(2Z + l)C^(co S7 ), (46) 



where 7 is the angle between the directions ei and e 2 and 
Pi is a Legendre polynomial. The coefficients Ci are the 
multipole moments. In what follows, we will be mainly 
interested in the so-called band power STi defined by the 
following expression 



8^ = T d 1(1 + 1) 



Ci 

Z7T 



(47) 



where Tq ~ 2.7 K. The band power has now been mea- 
sured on a wide range of angular scales from 10' to 90° 
corresponding roughly to I £ [2,700]. Almost 80 data 
points have been measured. Recently new data obtained 
by the balloon-borne experiments BOOMERanG 
and MAXIMA- 1 |§|f] have been published. They 
clearly show a detection of the first Doppler peak at the 
expected angular scale ~ 1° corresponding to the size of 
the Hubble radius at recombination. 

On the theoretical side, the multipoles moments de- 
pend on the initial spectra for scalar and tensor modes 
and on how the perturbations evolve from the initial time 
(after inflation) until now. This evolution is determined 
by the values of the cosmological parameters, i.e. by the 
value of the Hubble constant (h), of the total amount of 
matter present in our Universe (^o), of the cosmological 
constant (O^), of the baryons density parameter (f^b) 
and of the cold dark matter density parameter (f2 c dm)- 
Constraints already exist on some of these parameters. In 
particular, as already mentioned above, J7a — 0.7 accord- 
ing to the SNIa measurements and h 2 Q.h ~ 0.019 ± 0.002 
according to BBN f34| , [55|| . We also assume fi = 1 in 
agreement with the inflation paradigm which has been 
confirmed by the recent CMB anisotropy measurements. 
For the initial spectra, it is traditional to assume that 
they are of the power-law form 



k 3 P 9 (k) = A s k ns -\ k 3 P h (k) = A T k ri 



(48) 



where the scalar and tensor spectral indices ng and nx are 
related by ng — 1 = «t- This last equation is also valid for 
zeroth order slow-roll inflation. It should be noticed that, 
a priori, this choice is not the most relevant one since 
slow-roll inflation is certainly more physically motivated. 
For spectral indices close to ng = 1, we expect a small 
difference. This is no longer true for larger tilts. Inflation 
predicts the presence of gravitational perturbations and 
the tensor to scalar amplitude ratio is given by 



200 
~~9 



-7lT' 



(49) 



This equation is valid for power-law inflation with nx not 



too large^] or for zeroth order slow-roll inflation. A last 
remark is in order at this point. All the plots displayed in 
this article are COBE normalized in the following way: 
the position of the Sachs- Wolfe plateau is tuned such 
that it best fits the COBE data points. In practice, this 
almost amounts to normalize the spectrum to Cio- 

In this section, we first study the general properties 
of the multipoles moments in the quintessence cold dark 
matter model (QCDM) and point out the main differ- 
ences with the standard cold dark matter (sCDM) and 
the cosmic concordance model ( ACDM) . We also display 
the corresponding baryonic matter power spectra, given 
by 



\S(k)\ 2 ^ 



5pb 



Pb 



(50) 



which is the square of the Fourier transform of the 
baryonic density contrast. Then, we compare the pre- 
dictions of the QCDM model for the Ratra-Peebles 
and SUGRA tracking potentials with the COBE Q, 
BOOMERanG MM, MAXIMA-1 MM and Saska- 



toon 1 56 1 data. We do not attempt to perform a detailed 
statistical analysis but we rather indicate roughly how 
the different models can fit the observational data. 

We now turn to simple considerations about the shape 
of the CMB spectrum. The corresponding band power for 
the Ratra-Peebles and SUGRA potentials are displayed 
in Figs. I and | for h = 0.5, Q h = 0.05, Q Q = 0.7, 
fi c dm = 1 — — f^Q, ns — 0.99 and the tensor contribu- 
tion neglected. The former set of cosmological parame- 
ters has been chosen just for the sake of illustration and 
discussion. For simplicity, we start with a comparison of 
the quintessence multipole moments with those obtained 
in the ACDM model with similar cosmological parame- 
ters. Firstly, since fl m = f2 cc j m +^b is the same in the two 
models, the redshift of equivalence between matter and 
radiation z oq = f2 m /f2 r , where J7 r = f2 7 + Q^, is also the 
same in both cases. Therefore, the first peak is boosted in 
the same way by the early integrated Sachs- Wolfe effect 
(due to the time variation of the two Bardeen potentials 
during recombination, see |57j| ) and, a priori, one expects 
the same first peak height. Secondly, the dark energy 
component (cosmological constant or quintessence) has 
a negligible contribution before recombination and, as a 
consequence, the evolution of the perturbations before 
the last scattering surface is the same in the two mod- 
els (see the previous section). Thus, one expects again 
identical acoustic peak patterns. However, despite the 
previous considerations, the position of the peaks differs 
because the angular distance-redshift relation is modified 
at small redshift since the equation of state of the cos- 
mological constant and of quintessence is not the same. 



*For power-law inflation, the exact expression is given by 
A T /A s = -(200/9)n T /(l - n T /2). 
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FIG. 4. Multipole moments for the Ratra-Peebles potential 
for two values of a, a = 6 (red curve) and a = 11 (green 
curve) and with cosmological parameters equal to h = 0.5, 
n h = 0.05, D.q = 0.7, fi c d m = 1 - fib - «q, n s = 0.99, 
At = 0. The curves are compared with those obtained in 
the SCDM model (blue curve) and in the ACDM model (pink 
curve) . 
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FIG. 5. Multipole moments for the SUGRA potential for 
two values of a, a = 6 (red curve) and a = 11 (green 
curve) and with cosmological parameters equal to h = 0.5, 
fib = 0.05, D.q = 0.7, fi c d m = 1 - fib - fiq, n S = 0.99, 
At = 0. The curves are compared with those obtained in 
the SCDM model (blue curve) and in the ACDM model (pink 



The closest to —1 the equation of state parameter is, the 
largest the shift of the peaks to small angular scales is. As 
a consequence, the peaks in the ACDM model are more 
shifted to the right than in the QCDM model. Another 
feature is that the height of the first peak is not the same 
in the two types of scenarios. Indeed, at small redshift, 
the gravitational potential does not behave exactly in the 
same way in the two models especially because there are 
scalar field perturbations in the QCDM scenario. This 
results in a different contribution of the late integrated 
Sach- Wolfe effect j5?j] which affects the overall normal- 
ization of the spectrum. As a consequence, the height 
of the first peak is lower in the model which produces 
a strong late integrated Sachs-Wolfe effect, i.e. in the 
QCDM model. 

The exact shape of the quintessence potential also 
matters and different potentials lead to different CMB 
anisotropics. The SUGRA potential and the cosmologi- 
cal constant lead to very similar CMB anisotropy spec- 
tra, whereas the difference is stronger in the case of the 
Ratra-Peebles potential. This is mainly due to the fact 
that the equation of state parameter is generically closer 
to —1 in the first case than in the second one. Another 
difference is that the Ratra-Peebles potential produces a 
larger late integrated Sachs- Wolfe contribution than the 
SUGRA potential. This results in a different normal- 
ization for both models (note that the normalisation de- 
pends on a) which has for consequence different height 
of the first Doppler peak. Of course, this difference is 
also visible in the power spectrum at large scales. Maybe 
the most interesting property is the following one. The 
cosmic equation of state (almost) does not depend on a 
in the case of the SUGRA potential. Then, in the same 
manner, the CMB anisotropies do not depend on a con- 
trary to the case of the Ratra-Peebles potential. This 
means that the multipole moments displayed in Fig. ^ 
are a generic predictions of the SUGRA QCDM model. 

For the sake of completness, let us now describe the 
corresponding matter power spectra. They are displayed 
in Fig. ^| and The matter power spectrum also 

depends on the nature of the dark energy component 
(cosmological constant or quintessence) but the differ- 
ence between the cosmological constant scenario and 
a quintessence scenario is less important. The matter 
power spectrum shows a peak the location of which is 
given by the Hubble radius at equivalence. In the ACDM 
and QCDM scenarios, the peak is at the same location 
contrary to the sCDM case for which the peak is located 
at smaller scales. Also, in models with low matter con- 
tent, the ratio f2b/^cdm is higher which results in the 
presence of smooth oscillations at small scales. As for 
the CMB anisotropy spectrum, the small scales are sim- 
ilar in the ACDM and QCDM scenarios and important 
differences only occur on larger scales which are more 
affected by the change in the cosmic equation of state. 

Let us now study in more details and for more real- 
istic values of the cosmological parameters, the position 
and the height of the first Doppler peak. We start with 
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FIG. 6. Matter power spectrum for the Ratra-Peebles po- 
tential for two values of a, a — 6 (red curve) and a = 11 
(green curve) and with cosmological parameters equal to 
h = 0.5, fib = 0.05, = 0.7, Qcdm = l-fib-0 Q , n s = 0.99, 
^4t = 0. The curves are compared with those obtained in the 
SCDM model (blue curve) and in the ACDM model (pink 
curve) . 
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FIG. 7. Matter power spectrum for the SUGRA potential 
for two values of a, a = 6 (red curve) and a = 11 (green 
curve) and with cosmological parameters equal to h = 0.5, 
n h = 0.05, D.q = 0.7, fi c d m = 1 - n b - fiq, ns = 0.99, 
At = 0. The curves are compared with those obtained in 
the SCDM model (blue curve) and in the ACDM model (pink 
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the location of the first peak (denoted in what follows 
by li) and we study it in the plane (f2 m ,/i) with the 
following values of the other cosmological parameters: 
/i 2 £lb = 0.019 (the value predicted by standard BBN), 
Q A Q = 0.7, = and n s = 0.99. The case of the 
cosmological constant is displayed in Fig. |, the case of 
the Ratra-Peebles QCDM model in Fig. |9 and the case 
of the SUGRA QCDM model in Fig. |o[ These plots 
confirm the qualitative predictions made above and in 
particular the fact that, in general, > lf VGKA > lf F . 
If one assumes that f2 m ~ 0.3 (since we have assumed 
J7a — 0.7) and h ~ 0.62, this last value being consis- 
tent with the Hubble Space Telescope (HST) and SNIa 
measurements, then we obtain if ~ 225, Z^ugra ^ 2 20 
and if^ ~ 200. It is interesting to compare these values 
with the recent measurements of the first peak performed 
by BOOMERanG and MAXIMA- 1. The BOOMERanG 
data indicate that l\ = 197 ± 6 |f20| , ^lf which is compati- 
ble with the Ratra-Peebles potential and a spatially flat 
Universe. On the other hand, the MAXIMA-1 data are 
consistent with a first peak located at l\ ~ 220 |p2 23 
which is, this time, in agreement with a cosmological con- 
stant or the SUGRA QCDM model. 

Let us now study the height of the first Doppler peak. 
We study its variation in the plane (fibj^s) for the fol- 
lowing values of the cosmological parameters: h = 0.62, 
J7a,q = 0.7. The case of the ACDM model is displayed 
in Fig. [ll] whereas the cases of the Ratra-Peebles QCDM 
and SUGRA QCDM are presented in Figs. [l| and [l|, 
respectively. We would like to emphasize that the im- 
portance of gravitational waves is crucial in this case. In- 
deed, as already mentioned, the presence of gravitational 
waves modifies the normalization and, as a consequence, 
the height of the peaks. The BOOMERanG data indicate 
that 5T 20 o ~ 69 ± 8 /iK MM whereas the MAXIMA-1 
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Matter density £2 m 

FIG. 9. Contour plots of the first Doppler peak location 
in the (fl m ,/i) plane for the Ratra-Peebles QCDM case. The 
other cosmological parameters are h?Q.h = 0.019, Slq = 0.7, 
A T = and n s = 0.99. 



Baryon density Cl b 

FIG. 11. Contour plots of the height of the first peak in 
the (flb, ns) plane with h = 0.62, = 0.7 for the case of the 
ACDM model. 
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FIG. 10. Contour plots of the first Doppler peak location 
in the (fi m , h) plane for the SUGRA QCDM case. The other 
cosmological parameters are h 2 Q h = 0.019, Qq = 0.7, At = 
and n s = 0.99. 
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FIG. 12. Contour plots of the height of the first peak in 
the (f2b,ns) plane with h — 0.62, S1q = 0.7 for the case of the 
Ratra-Peebles QCDM model. 
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FIG. 13. Contour plots of the height of the first peak in 
the (fib) ns) plane with h = 0.62, S7q = 0.7 for the case of the 
SUGRA QCDM model. 
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FIG. 14. Contour plots of the location of the second peak 
in the (Cl m ,h) plane with h = 0.62, Q A = 0.7, h 2 Q b = 0.019, 
ns = 0.99 for the cosmological constant case. 



ones give ST 220 — 78 ± 6 [iK [£2|,g3| , this discrepancy be- 
ing possibly explained by problems in the calibration of 
these experiments. If we adopt the value fib — 0.0595, 
compatible with BBN, we see that, in the Ratra-Peebles 
QCDM model, a height of the first peak compatible with 
the BOOMERanG and MAXIMA- 1 data leads to a value 
of the scalar spectral index such that rig > 1. This is not 
compatible with standard inflation and cannot be real- 
ized with one scalar field. We interpret this as a new 
evidence that the Ratra-Peebles QCDM model (at least 
with this value of a) is excluded by the observations. For 
the cases of ACDM and SUGRA QCDM, we learn from 
the previous plots that the spectral index must be very 
close to one. 

We now turn to the study of the second Doppler peak. 
First of all, we should say something about the obser- 
vational situation. With regards to the detection of a 
second peak, it is difficult to deduce something from the 
BOOMERanG data. The error bars are still large and the 
data are, for the moment, compatible with a second peak 
(with a height maybe smaller than predicted by standard 
inflation) but also with no peak at all, even if one can see 
a small rise of the signal at l 2 550 Only 5% of 

the data of this experimement have been analysed so far 
and one should wait for the rest of the data analysis to 
be completed. On the other hand, the MAXIMA-1 show 
"a suggestion of a peak at l 2 — 525" (22) the height of 
which would be ST525 — 48 f^K. One could even argue 
that the beginning of a third peak has been observed. In 
fact, considering all the uncertainties in such measure- 
ments, we are of the opinion that a reasonable attitute 
is simply to wait for more data. On the theoretical side, 
it was argued by Kamionkowski and Buchalter |l0| that 
the location of the second peak can probe the dark en- 
ergy density. The main idea is to study the contour plots 
of I2 in the plane (f2 m ,/i). Then, a measurement of I2, 



knowing h by other means , immediately determines the 
value of f£ m . It was claimed in Ref. Q that this strategy 
does not depend on whether the dark energy is a cosmo- 
logical constant or a quintessence field. We show that 
this claim is not correct and that the nature of the dark 
energy matters. The contour plots of l 2 in the case of 
a cosmological constant are displayed in Fig. [li] for the 
cosmological parameters given by h = 0.62, fljy = 0.7, 
/i 2 f2b = 0.019, ns = 0.99. These plots are in agreement 
with the results found in Ref. j4(|. The corresponding 
contour plots for the Ratra-Peebles and SUGRA CDM 
models are presented in Figs. 15 and [l^. In addition, in 
order to show that there is indeed an important differ- 
ence, we also display the contour plots for a cosmological 
constant which, for a given value of I2, is always above 
the QCDM curve. The fact that there is a difference 
does not totally invalidate the idea of Ref. Q . But it 
means that, in order to use it, we should first identify 
the physical nature of the dark energy, for example with 
a measurement of its equation of state parameter. As 
for the first peak, we have l 2 > /f UGRA > lf F . Roughly 
speaking, for h ~ 0.62, fl m ~ 0.3, we have l 2 — 550, 
/SUGRA ^ 525 an( j ;RP ^ 500. Interestingly enough, 
the SUGRA QCDM model seems to predict the correct 
location of the "suggested second peak" ^2|, just in be- 
tween the location predicted by the ACDM model and 
the Ratra-Peebles QCDM model. Of course, it is prema- 
ture to conclude and only more data could allow to know 
whether this is indeed the case or whether this is just a 
coincidence. 

Finally, we display the multipole moments for the 
ACDM model, the Ratra-Peebles QCDM model and the 
SUGRA QCDM model in Figs. 0, [l| and ||, respec- 
tively, for the following cosmological parameters (de- 
duced from the previous considerations): h — 0.62, 
tt A = 0.7, tt h = 0.0595 and n s = 0.99. The data points of 
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FIG. 15. Contour plots of the location of the second peak 
in the (fi m , h) plane with h = 0.62, fi A , Q = 0.7, h 2 Q b = 0.019, 
n s = 0.99 for the Ratra-Peebles QCDM model. The corre- 
sponding contour plots for the cosmological constant (upper 
curves) are also displayed for comparison. 
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FIG. 17. Band power ST t for the ACDM model with 
h = 0.62, D.a = 0.7, fib = 0.595 and n s = 0.99. The data 
point are those of the COBE, BOOMERanG, MAXIMA-1 
and Saskatoon experiments. 
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FIG. 16. Contour plots of the location of the second peak 
in the (fi m , h) plane with h = 0.62, Q. a ,q = 0.7, h 2 Q b = 0.019, 
n s = 0.99 for the SUGRA QCDM model. The corresponding 
contour plots for the cosmological constant (upper curves) are 
also displayed for comparison. 
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FIG. 18. Band power 5T t for the Ratra-Peebles QCDM 
model with h = 0.62, Q.q = 0.7, fi b = 0.0595 and n s = 0.99. 
The data point are those of the COBE, BOOMERanG, MAX- 
IMA-1 and Saskatoon experiments. 
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FIG. 19. Band power 8Ti for the SUGRA QCDM model 
with h = 0.62, D.Q = 0.7, D. h = 0.0595 and n s = 0.99. The 
data point are those of the COBE, BOOMERanG, MAX- 
IMA- 1 and Saskatoon experiments. 



COBE, BOOMERanG, MAXIMA-1 and Saskatoon have 
been added to the plots for comparison. These curves 
represent the predictions of each model and special at- 
tention must be paid to third peak which is certainly one 
of the next important experimental challenge. In Fig. |2^, 
we present the three curves together in order to make the 
comparison easier. It should be emphasized again that 
the multipole moments predicted by the SUGRA QCDM 
model are unique in the sense that they do not depend on 
the free parameter in the potential. From these plots, we 
see that the SUGRA QCDM model is, among the three 
models studied here, the best fit of the MAXIMA-1 data. 
It is the only model for which the theoretical curve 5T 
versus I goes through all the la error bars of this experi- 
ment. However, we should be careful not to overestimate 
the relevance of this result since uncertainties are still 
large, for instance because the comparison of the calibra- 
tions of different experiments is always a difficult task. 
We should also keep in mind that 2a deviations are al- 
ways possible. Thus, we are waiting eagerly for the new 
data to see whether quintessence, and especially SUGRA 
quintessence, can confirm the hints of this article and fits 
the data better than the other QCDM models. 



V. CONCLUSION 
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FIG. 20. Band power 8Ti for the ACDM model (blue 
curve), Ratra-Peebles QCDM model (green curve) and the 
SUGRA QCDM model (red curve) with h = 0.62, Q A ,Q = 0.7, 
fib = 0.0595 and ns = 0.99. The data point are those of the 
COBE, BOOMERanG, MAXIMA-1 and Saskatoon experi- 
ments. 



The quintessence scenario provides a general frame- 
work within which the issue of the energy density of the 
Universe can be tackled. In particular long-standing is- 
sues such as coincidence problem (and maybe the fine- 
tuning problem) receive reasonnable answers for a class 
of models possessing the property of tracking fields, i.e. 
the evolution of the quintessence field is driven at small 
redshift towards an attractor independently of the ini- 
tial conditions. In the same spirit it seems very entic- 
ing to draw consequences of the quintessence hypothesis 
on other cosmological observables, the most prominent 
ones being the cosmological anisotropies. Recent mea- 
surement of the CMB anisotropies by the BOOMERanG 
and MAXIMA-1 experiments give a first indication on 
the location of the peaks in the CMB multipoles. It 
seems therefore topical to understand the consequence 
of the quintessence hypothesis on the CMB anisotropies. 

In this paper we have confronted analytical methods 
with numerical results. Using the former we establish 
that the quintessence perturbation are independent of 
the initial conditions. This is confirmed by a full nu- 
merical computation. This allows us to study the CMB 
anisotropies. In particular we have paid particular at- 
tention to the comparison between three possible mod- 
els: the cosmological constant model, the Ratra-Peebles 
and SUGRA quintessence models. We have also com- 
pared these three models with the existing data from the 
BOOMERanG and MAXIMA-1 experiments. As a rule 
the location of the first peak is shifted to the right for 
models having an equation of state u> closer to —1. This 
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entails that the location of the first peak for the first peak 
of the MAXIMA- 1 data is fitted by the SUGRA model. 
Similarly the location of the second peak around I2 — 525 
as suggested by MAXIMA- 1 seem to indicate that the 
SUGRA model comes closer to be the best of these three 
models. One of the foreseeable challenges will be to carry 
out a thorough analysis of the forthcoming data in order 
to distinguish these three models even more clearly. 

From the particle physics point of view most of the 
quintessence models discussed so far have neglected the 
crucial effects of SUSY breaking. It may well be that 
the effects of SUSY breaking, on top of necessitating a 
severe fine-tuning of the cosmological constant, will in- 
duce drastic modifications in the functional form of the 
quintessence potential. It is certainly a tantalizing chal- 
lenge to include the effects of SUSY breaking within the 
supergravity models of quintessence . On the other 
hand there exists the possibility that the cosmological 
constant problem will be resolved using ideas stemming 
from extra-dimension scenarios involving an effective su- 
persymmetry in four dimensions . The investigation of 
such models might well shed new light on the origin of 
the quintessence field. 

As must be clear by now the issues raised by the cos- 
mological constant problem, the quintessence scenario 
and its proper understanding within particle physics are 
many-fold. The experimental results which will be avail- 
able in the near future might help in disentangling some 
of these very conspicuous matters. 
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